arXiv:1505.08103v2 [physics.optics] 17Jun2015 


Selective coupling of optical energy into the fundamental diffusion mode of a 

scattering medium 

Oluwafemi S. Ojambati, Hasan Yilmaz, Ad Lagendijk, Allard P. Mosk, and Willem L. Vos 
Complex Photonic Systems (COPS), MESA-h Institute for Nanotechnology, 

University of Twente, P.O. Box 217, 7500 AE Enschede, The Netherlands 
(Dated: June 17th, 2015) 

We demonstrate experimentally that optical wavefront shaping selectively couples light into the 
fundamental diffusion mode of a scattering medium. The total energy density inside a scattering 
medium of zinc oxide (ZnO) nanoparticles was probed by measuring the emitted huorescent power 
of spheres that were randomly positioned inside the medium. The huorescent power of an optimized 
incident wavefront is observed to be enhanced compared to a non-optimized incident wavefront. 

The observed enhancement increases with sample thickness. Based on diffusion theory, we derive 
a model wherein the distribution of energy density of wavefront-shaped light is described by the 
fundamental diffusion mode. The agreement between our model and the data is striking, not in the 
least since there are no adjustable parameters. Enhanced total energy density is crucial to increase 
the efficiency of white LEDs, solar cells, and of random lasers, as well as to realize controlled 
illumination in biomedical optics. 


Numerous physical phenomena are described by dif¬ 
fusion [l|-0. Diffusion is a process that leads to uni¬ 
form spreading of matter or energy as a result of ran¬ 
domness ii. Diffusion theory accurately describes the 
propagation of the energy density of multiply scattered 
waves in disordered scattering media (^. Isl. [iq4i^ . Upon 
averaging over the disorder, waves become diffuse after 
a distance of the order of one transport mean free path 
£ and the energy density of the waves acquires a typical 
shape, shown in Fig. [Hb). The derivative of the energy 
density at the exit surface is related to the transport of 
energy, and yields the waves-equivalent of the well-known 
Ohm’s law, T ^ £)L, where L is thickness of the scatter¬ 
ing medium. 


In a slab geometry, the solution of the diffusion equa¬ 
tion can be expressed as a sum over a complete set 
of eigensolutions with imaginary frequency [l^- In 
Fig. UTa), we show the first three eigensolutions. When 
a plane wave is incident on a scattering medium, energy 
is coupled into all eigensolutions, which sums up to give 
the non-optimized energy density Wd shown in Fig.[T]^b). 
A fundamental question we seek to address is the oppor¬ 
tunity of changing the internal energy by selectively cou¬ 
pling energy only into the fundamental diffusion eigen¬ 
mode with index m = 1 shown in Fig. HJa). It is of par¬ 
ticular interest when the total energy coupled into the 
fundamental diffusion mode is greater that of the non- 
optimized energy density Wd as shown in Fig. In 

some particular cases, the fundamental diffusion eigen¬ 
mode has a greater total energy density than the un¬ 
optimized energy density Wd as shown in Fig. [IJb). In 
the case of light, which is the subject of our work, an 
enhanced energy density inside the scattering medium 
is important for applications, such as enhanced energy 
conversion in white LEDs [iMi, efficient light harvest- 
ing in solar cells [2 qI-[ 2^. low-threshold random lasers 


23|, 12^, and controlled illumination in biomedical 



FIG. 1. (Color) (a) The first three eigensolutions of the dif¬ 
fusion equation, where m is index of the eigensolution. (b) 
The energy density of optimized light and non-optimized light 
are shown as the red and blue curves respectively. The en¬ 
ergy density is reduced with the diffusion constant D, the 
incident intensity E and the effective thickness of the sample 
Lex = L + Zei + ^e 2 , ^ei and Zq 2 are the extrapolation lengths 
at the front and back surfaces of the sample respectively. 


optics |25j . 

The total transmitted intensity through a scattering 
medium can be made to differ from Ohm’s law by wave- 
front shai 
jugation [T 
trix 


26|-[3l| , time reversal [32|, [33| , phase con- 
I , and control based on transmission ma- 
361. l37j . In wavefront shaping, the spatial phase of 


the incident field on the scattering medium is controlled 
in order to enhance the intensity in a diffraction-limited 
spot at the back surface of the sample. Only numerical 
calculations 38I44Q| and a single-realization experiment 
of elastic waves with a shaped incident wavefront 
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have been used to study the distribution of energy den¬ 
sity inside a two-dimensional (2D) scattering medium. 
The distribution observed in these calculations and in 
the single-realization experiment is a symmetric func¬ 
tion peaked at the middle of the sample, which is sim¬ 
ilar to the fundamental eigenmode m = 1. The change 
in the energy density has so far not been experimen¬ 
tally observed inside a three-dimensional (3D) scattering 
medium. 

In this Letter, we demonstrate experimentally the se¬ 
lective coupling of light into the fundamental mode of the 
diffusion equation by using wavefront shaping. We probe 
the total internal energy, which is integral of the position- 
dependent energy density inside a 3D scattering medium. 
As a probe, we employ fluorescent spheres randomly po¬ 
sitioned inside the medium. We observe that the total 
energy increases when the incident light is shaped. The 
enhancement in fluorescent power increases with sample 
thickness. To interpret our results, we propose a model 
wherein the energy density of wavefront-shaped light is 
described by the fundamental eigensolution of the diffu¬ 
sion equation. Our model has no adjustable parameters 
and agrees well with the experimental results. 

In our experiments, we study a scattering medium, 
which is a layer of spray-painted zinc oxide (ZnO) 
nanoparticles on a microscope glass slide of thickness 
0.17 mm. The transport mean free path of similar sam¬ 
ples was determined from total transmission measure¬ 
ments to be 0.6 ± 0.2 //m [42|. Inside the ZnO samples, 
dye-doped polystyrene spheres with diameter 50 nm are 
randomly dispersed, as illustrated in Fig. O The fluo¬ 
rescent spheres are excited by incident laser light with 
wavelength Ai = 561 nm and emit fluorescent light at a 
different wavelength A 2 = 612 nm. In order to ensure 
that the spatial distribution of the energy density at Ai 
inside the scattering medium is not perturbed by the ab¬ 
sorption from the probing fluorescent spheres , we use 
samples with a low density of spheres and with a high 
albedo ji^ . 

A phase-only liquid crystal spatial light modulator 
(SLM) (Holoeye Pluto) shapes the wavefront of the laser 
light incident on the sample, such that the intensity is 
focused in a diffraction-limited spot at the back surface 
of the sample. We used the piece-wise sequential algo¬ 
rithm described in Ref. to find an optimized incident 
wavefront. The back surface of the sample is imaged to 
the chip of an electron multiplying charged-coupled de¬ 
vice (EMCCD) camera to collect the total fluorescent 
intensity at A 2 . A combination of a dichroic mirror, a 
low-pass filter and a notch filter blocks the incident light 
at Ai from reaching the EMCCD ji^ . 

To obtain ensemble-averaged data that can be com¬ 
pared to theory, we need to average over different realiza¬ 
tions of scatterers in the sample. We therefore performed 
automated sequences of wavefront shaping measurements 
while between two consecutive measurements, the sample 
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FIG. 2. (Color) Schematic drawing of the method to probe 
the total energy density inside a scattering medium. The 
scattering medium is an ensemble of disordered ZnO particles 
in air. The medium is illuminated with a shaped incident 
wavefront such that the incident light at Ai (green intensity) 
is optimized on a diffraction-limited spot at the back of the 
sample. The scattering medium is lightly doped with fluores¬ 
cent spheres randomly positioned inside the medium to probe 
the energy density inside the sample. The total fluorescent 
power emitted from the fluorescent spheres at A 2 (red inten¬ 
sity) is measured by EMCCD. 


was translated to a different realization by a piezo stage. 
In each measurement, we measured the total fluorescent 
power with the optimized pattern on the SLM and 
then measured the total fluorescent power with an 
incident wavefront optimized for a different uncorrelated 
position. We define the fluorescent power enhancement 
r]f as the ratio of the two fluorescent powers, = P^/P^. 

We determine the fidelity jqp that quantifies the over¬ 
lap of the experimentally generated field with the ideal 
controlled field [ 2 ^. The fidelity jyp achieved in a spe¬ 
cific experimental run can be obtained by dividing the 
optimized power in the target spot by the average total 
transmitted power without optimization. We performed 
100 wavefront shaping experiments, each on a different 
position on a L = 22.8 /im ± 0.95 thick sample. In each 
experiment, we determine the fidelity jyp and fluorescent 
power enhancement r]f . Eactors such as inhomogeneity of 
the sample thickness, measurement noise and instability 
in environmental conditions result into variation of the 
fidelity I 7 P [H, [ 4 ^. Although these factors are undesir¬ 
able, they have the advantage of giving a wide range of 
I 7 P to investigate. 

In Eig.[3]we show the enhancement in fluorescent power 
versus the fidelity jyp. Interestingly, we see an en¬ 
hancement in the fluorescent power by up to about 10 % 
as I 7 P increases to about 0.035. This increase implies 
that the total energy density for optimized incident wave- 
fronts is higher than the total energy density of unopti¬ 
mized incident wavefronts. If the spatial distribution of 
energy density would be unmodified by wavefront shap¬ 
ing, then the fluorescent intensity enhancement would 
have been constant at 1 , which is obviously not the case. 
The slope of the linear regression fit to the data is 3.6 
with a standard error of 0.2, and upper and lower 95% 
confidence intervals of 4.1 and 3.2, respectively. The cor- 
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FIG. 3. (Color online) Measured fluorescent power enhance¬ 
ment ? 7 f versus the fidelity for an L — 22.8 /xm zb 0.95 
thick ZnO sample. The black squares are 100 experimental 
data points obtained at different positions on the sample. The 
solid red curve is a linear regression through the data and the 
blue dashed curves are the 95% confidence interval. The green 
dash-dotted curve is the expected curve if the distribution of 
light with optimized wavefronts were the same as with diffuse 
light. 

relation coefficient r of the data is 0.9, which confirms 
that our data show a linear trend. The measured fluores¬ 
cent power enhancement rji has contributions from both 
the perfectly shaped wavefront and from the background 
intensity, which is the uncontrolled part of the intensity. 
We therefore express r]i in terms of the fidelity jyp as 
rji =7^f|7p + (l — | 7 p), where is the fluorescent power 
enhancement extrapolated to the limit of perfect fidelity 
| 7 p ^ 1. The second term is the contribution from the 
background intensity. For the result shown in Fig. O we 
find = 4.6 ± 0.48. 

We studied samples with thicknesses L ranging from 2 
/im to 22 jam and on each sample we performed 100 to 
130 wavefront shaping experiments. Since the fidelity jyp 
decreases with increasing sample thickness [HI, we de¬ 
rived for each sample the extrapolated fluorescent power 
enhancement rj^ to allow for a comparison between sam¬ 
ples. In Fig.m we show that the extrapolated fluorescent 
power enhancement r]f increases with sample thickness 
L, which means that wavefront shaping serves to opti¬ 
mally store energy in the volume of the medium. The un¬ 
certainty in the extrapolated fluorescent power enhance¬ 
ment increases with sample thickness, since the fidelity 
decreases for thicker samples. The horizontal error bars 
denote the standard deviation of the measurement of the 
sample thickness on different positions on the sample. 
For perfect fidelity, the total fluorescent power inside a 
22.8 jam thick sample is 4.6 times greater than the total 
fluorescent power for non-optimized light. 

To interpret our experimental results, we employ dif¬ 
fusion theory [dj. We obtained the diffuse energy den¬ 
sity Wd shown in Fig. [TJb) from the diffusion equation. 
For light with an optimized incident wavefront, the dis¬ 
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FIG. 4. (Golor online) Fluorescent power enhancement r]f 
in ZnO scattering samples versus sample thickness. The red 
squares are the measured fluorescent power enhancement ex¬ 
trapolated to unity intensity control. The blue solid line is the 
calculated fluorescent power enhancement from Eq. [U The 
green dash-dotted curve is for an invariant distribution of en¬ 
ergy density along the sample depth. 


tribution of the energy density Wq inside the medium 
is a-priori unknown. With the optimized phase inci¬ 
dent on the sample, light is coupled to the transmis¬ 
sion eigenmodes of the wave equation with the highest 
transmission. Since both the wave equation and the dif¬ 
fusion equation describe the same physical system, we 
expect that the ensemble-averaged energy density of the 
transmission eigenmodes with the highest tranmission is 
equivalent to the diffusion eigensolution that contributes 
the most to the total transmission. We show in Fig. [5] 
the contribution to the total transmission [HI of the first 
six eigensolutions. The fundamental eigensolution m = 1 
contributes the most to the total transmission, even more 
than the total transmission, which is a summation of con¬ 
tribution of all the eigensolutions. We therefore hypoth¬ 
esize that the energy density distribution of optimized 
light is identical to the fundamental eigensolution of the 
diffusion equation. The validity of this hypothesis is ver¬ 
ified when we compare our model to experimental data. 

It has been shown experimentally in Ref. [Hj and the¬ 
oretically in Ref. [49| that the total transmission of opti¬ 
mized light is equal to Tq = 2/3. We therefore scale the 
energy density of wavefront-shaped light such that the 
total transmission is equal to Tq = 2/3. In Fig. [I{b), we 
show the scaled energy density of light with an optimized 
wavefront that clearly deviates from the distribution of 
diffuse light. In addition, Davy et al. theoretically cal¬ 
culated the internal energy density distribution of trans¬ 
mission channels with a transmission coefficient of unity 
and found a parabolic solution [HI- Since the boundary 
conditions used in Ref. [HI only apply to a medium that 
is index-matched to the surrounding media, the model 
does not pertain to our experiments. Nevertheless, it is 
interesting to note that the parabolic function found by 
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Davy et al. and our sine function are both symmetric 
functions peaked in the middle of the sample. 


FIG. 5. (Color online) The contribution to the total trans¬ 
mission of the first six eigensolutions are represented by the 
red bars while the total transmission, which is the sum of all 
transmissions of the individual eigensolutions is represented 
by the blue bar. 


Taking into account the diffusion of the emitted fluo¬ 
rescent light propagating through the sample we analyt¬ 
ically model the fluorescence enhancement as 


v!{L) = 


2ie"xsec(^) TT^ei cos(^) - rex[sin(^) - sin(f^)] 


cos(^) 


Stt^ 


Lz[^.[L^-\-3L(Zel-\-Ze2)+6ZelZe2] 

61 /ex 




- z' 

Gnj^inj 


( 1 ) 


where Lex = L-\-Zei+Ze 2 is the effective sample thickness, 
L' = L + ^^ei, and z[^- = Zini + Zei^ and zinj is the injection 
depth at which the incident light becomes diffuse and 
it accounts for the angular distribution of the incident 
shaped wavefront [H^ . In order to compare our model to 
our experimental results, we plot in Fig. |4] the analytic 
model for rjf versus sample thickness L. Our model agrees 
well with our experimental result. There are no freely ad¬ 
justable parameters in our model. If the spatial distribu¬ 
tion of both wavefront-shaped and unwavefront-shaped 
light would have been the same, then r]f would been con¬ 
stant equal to 1 for all sample thicknesses as shown in 
Fig.m which does not agree at all with our observations. 
The excellent agreement between our model and our ex¬ 
perimental results confirms the validity of our hypothesis 
that the distribution of wavefront-shaped light inside the 
medium is modified, and that energy has been coupled 
into the fundamental diffusion mode. 

In our experiments we obtain the fluorescent power 
enhancement 77 / rather than the energy density enhance¬ 
ment 77ed- Therefore we define the enhancement of the 
energy density to be r]e = where and are 

the energy densities for optimized light and unoptimized 
light, respectively, both integrated over the whole sample 
thickness. We obtain 

r?ed(L) = ^r/f(L) + O(i/0, (2) 

where 0{L/l) includes higher orders of the series expan¬ 
sion in terms of L//, see supplementary material Q. We 


see from Eq. [2] that the total fluorescent power depends 
on the total energy density inside the medium. There¬ 
fore the observed increase of the fluorescence is indeed a 
measure of the increase of the energy density. 

In summary, we have experimentally demonstrated 
and theoretically modeled the selective coupling of light 
into the fundamental diffusion mode, by increasing the 
total energy density inside a scattering medium of ZnO 
nanoparticles by using wavefront shaping. Our results 
apply to other wave control methods in scattering me¬ 
dia, such as time reversal, phase conjugation, and control 
based on transmission matrix as well as to other types of 
waves such as microwaves, acoustic waves, elastic waves, 
surface waves, and electron waves. We expect our results 
to be relevant for applications that require enhanced to¬ 
tal optical energy density such as efficient light harvesting 
in solar cells especially in near infrared where silicon has 
low absorption; for enhanced energy conversion in white 
LEDs, which serves to reduce the quantity of expensive 
phosphor; for low threshold and higher output yield of 
random lasers; as well as in homogeneous excitation of 
probes in biological tissues. Last but not least, it will 
be fruitful to investigate possible relationships between 
the fundamental diffusion eigensolution and the univer¬ 
sal diffusion time obtained in Refs, d, IHlI. 1^. 
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